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Abstract 

A systematic analysis of transverse momentum distribution of hadrons 
r*| . produced in ultra-relativistic p + p collisions is presented. We inves- 

tigate the effective temperature and the entropic parameter from the 
non-extensive thermodynamic theory of strong interaction. We con- 
clude that the existence of a limiting effective temperature and of a 
limiting entropic parameter is in accordance with experimental data. 

1 Introduction 

General aspects of strong interactions up to center-of-mass energies \/s ~ 
10 GeV are well understood in terms of a self-consistent theory based on 
the Boltzmann-Gibbs statistics pQ. Hagedorn's theory establishes a connec- 
tion between the mass spectrum of highly excited hadrons and the density of 
states for fireballs, and provides correct descriptions for transverse momen- 
tum distributions and multiplicities of secondaries. 

For y/s > 10 GeV, however, theoretical and experimental results diverge. 
A generalized formalism was proposed taking into account the non-extensive 
thermodynamics [2J [3]. This generalization recovered the agreement between 
theory and experiment. 

Recently it was shown |3] that a self-consistent theory for hadronic sys- 
tems based on the non-extensive thermodynamics exists if, for x — > oo, 

p(x) -> 7 x- 5/2 e^ (1) 
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and 



a(x) -»■ bx a e f t oX , (2) 



where pirn) is the mass spectrum of hadrons, <j(E) is the density of states 
for hadronic systems, and a is given by 

- lV ° 0) 



2tt 2 /3 3/2 ' 

with 7 and b being constants and V D being the interaction volume. Here e x 
is the q-exponential function given by 

e* = [l + (g-l)a;]i^. (4) 

An important consequence of the self-consistency principle is the existence 
of a limiting effective temperature, T Q = l//3 , and of a limiting entropic 
parameter, q Q . 

In this work we perform a systematic analysis of experimental data on 
p+p collisions to verify if there exist T Q and q Q that allow a correct description 
of all experimental data. The set of experiments used in the present analysis 
is summarized in Table [TJ and covers a wide energy range. 



Table 1: Set of experimental data for p + p collisions. 



Experiment 


Energy (GeV) 




\v\ 


CMS (LHC) 


7000 [5] 


\v\ 


< 2,4 


CMS (LHC) 


2360 [5] 


\v\ 


< 2,4 


CMS (LHC) 


900 [5] 


\v\ 


< 2,4 


ALICE (LHC) 


900 [BJ 


\v\ 


< 0,8 


ATLAS (LHC) 


900 [7j 


\v\ 


< 2,5 



RHIC (BNL) 200 [8] 3, 3 < rj < 5, 



2 Transverse momentum distribution in non- 
extensive statistics 

A direct method to obtain the effective temperature of the hadronic system 
produced in hadron-hadron collisions is through the study of the transverse 
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momentum (jp?) distribution of secondaries. According to the non-extensive 
formalism proposed in Ref . [2] , the p^-distribution is given by 



1 da 



cp T 



l + (q-l)P\ pl+p 2 T + p 



-9/(9-1) 



dpi 



a dpx 

where pi is the longitudinal momentum and p is the mass of the hadron. 



(5) 



a- 



ALICE: 0.9TeVx 10 4 

RHIC: 0.2TeVx 10" 10 

ATLAS: 0.9TeVx 10~ 8 

CMS: 0.9TeVx 10~ 6 

CMS: 2.36 TeVx 10~ 2 

CMS: 7TeVx 10° 
Hagedorn 

Hagedorn-Tsallis 
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Figure 1: Fittings of Eq. [6] to experimental p^-distributions. 



An useful approximation was proposed in Ref. j3] assuming that pjph 
resulting 



1 da 



a dp< 



T 



i2(q - l)}- 1 '^ 



2q-l 



it 



3/2 



[1 + (q ~ 1H 



i 

9-1 



(6) 



where u = Ppt, and B(x,y) is the Beta-function. The dependence of the 
distribution on the entropic factor, q, and on the temperature, given by j3 = 
1/T, enable us to obtain both parameters by fitting the above expression to 
experimental data. This procedure has been already used by many authors 0, 
El El El [7J |8] . It is important to notice that Equation El should be regarded 
as an approximation to the distribution one would obtain from the Tsallis 
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entropy by using the thermodynamic relations [U E]> but to the use made 
here the approximation is good enough. 

Concerning the approximation on the mass, with the assumption that 
I^/Pl ~ 0, it should be mentioned that it can be avoided [9] with the restric- 
tion that one can not apply the formula obtained for general hadron (h + and 
h~), but only on specific particle distributions. Since in this work we aim to 
make a systematic analysis, and due to the larger number of information on 
(h + + /^-distributions, we opt to use the approximated formula El Even be- 
ing a good approximation, it causes the effective temperature obtained to be 
slightly shifted to higher values. This effect can be easily observed in Figure 
5 of Ref . [9] , where it can be seen that the peaks of the py-distributions for 
heavier particles are shifted to higher values with respect to that for pions. 

Applying the method described above to p+p collisions we get nice fittings 
to the experimental pr-distributions, as shown in Figure [TJ We observe that 
Eq. [6] correctly describes all data for pt up to 18 GeV/c. From these fittings 
we obtain the results for T and q that are shown in Fig [2j The temperature 
varies inside a relatively narrow range between 70 MeV and 90 MeV for all 
collisions with center-of-mass energy from 0.2 TeV up to 7 TeV. For energies 
above 0.9 TeV the temperature can be considered constant with T ~ 73 MeV. 
Also in the case of the parameter q the variations are relatively small in the 
energy range studied, and above 0.9 TeV it is approximately constant with 
q ~ 1.13. 

It is important to notice that T and q in Eq. E] are not completely in- 
dependent. There is a strong correlation between the two parameters, as 
can be observed in Fig. [3] where we show the 2D-plot of \ 2 distribution for 
the fittings shown in Fig. [TJ We clearly see the ellipses that evidences the 
correlation between the two parameters. Due to their correlation, as we vary 
T and q simultaneously along the line corresponding to the major axis of the 
ellipse, the \ 2 remains practically unchanged. This means that it is possible 
to obtain good x 2 with pairs (T, q) near the optimum point found in the 
fitting, indicated by crosses in Figure [3] 

A linear behaviour between T and q was already predicted by Wilk and 
Wlodarczyk [TO] [TT] , who proposed the relation 

T = T + (q-l)c, (7) 

where c is a constant depending on the energy transfer between the source 
and its surroundings and on thermodynamical properties of the medium [12] . 
Since T — > T as q — > 1, T is considered to be the Hagedorn's temperature. 

Assuming that Eq. [7] establishes a causal relation between T and q, then 
correlation is a necessary consequence [13]. Also, it is possible to obtain 
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Figure 2: The values for the parameters corresponding to the best fitted 
curve to pr-distributions. (a) T as a function of center-of-mass energy, (b) 
q as a function of center-of-mass energy. 



the regression coeficient, c, from the correlation. From the ellipses observed 
in Fig. [3 we obtain T =(192±15) MeV and c =-(0.95±0.10) GeV. It is 
interesting to notice that the value for T Q found here is in agreement with 
the critical temperature obtained in lattice QCD calculations [b£| IT5] . 

3 Testing the constant T hypothesis 

Now we slightly modify the fitting procedure by adopting a fixed value for 
T, as suggested by the self-consistency principle [1] and by the results shown 
in Fig. [2j We used different values for T from 60 MeV up to 120 MeV. Some 
results are shown in Fig. H] and we observe good fittings of Eq. [B] to the data 
in all cases, now with only q being a free parameter. 

These results show that we can fit all experimental data for p + p, for y^i 
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Covariance Matrix 



Covariance Matrix 



:.1S62 -0.001069 0.0005506 
001069 8.043e-06 -4.632e-06 



0.0005506 -4.662e-06 3.253e-0 



1.092 



X 2 = 5.460 
NDF (Degrees of freedom) = 30 
Z 2 ed = Z 2 /NDF = 0.182002 



Correlation = -0.9153 

1 - 4-, , . i , , , i , , , i , , , i , , 

0.076 0.078 0.08 0.082 0.084 O.C 



0.088 0.09 0.092 
T [QeV] 



cr1.125r 
1.1245: 

1.124 ; 
1.1235| 

1.123| 
1.1225| 

1.122; 
1.1215; 

1 .121 : 



0.023B3 -7.302e-05 5.695e-05 
-7.302e-05 3.459e-07 -3.77Be-07 



5.695e-05 -3.7786-07 5.29e-07 



X 2 = 96.983 
NDF (Degrees of freedom) = 44 
X 2 ed = r7NDF = 2.204159 
Correlation = -0.8833 



0.07 0.0705 0.071 0.0715 0.072 0.0725 0.073 

T [QeV] 



(a) RHIC: 200 GeV 



(b) ALICE: 900 GeV 



1.113 
1.112 
1.111 
1.11 
1.109 
1.108 
1.107 
1.106 
1.105 



0.3712 -0.001453 0.0007991 
-0.001453 6.149e-06 -3.759e-06 
0.0007991 -3.759e-06 2.639e-06 



X 2 = 8.016 
NDF (Degrees of freedom) = 31 
Z 2 ed = 2C 2 /NDF = 0.258585 
LCorrelation = -0.9333 



0.082 0.084 0.086 0.088 0.09 0.092 0.094 

T [GeV] 



Covariance Malrix 



0.2275 -0.0005596 0.0004172 
-0.0005596 2.399e-06 -2.4S4e-06 
0.0004172 -2.484e-06 2.997e-06 



X 2 = 36.931 
NDF (Degrees of freedom) : 
X 2 ed = rVNDF = 1 .846574 
Correlation = -0.9264 
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(c) ATLAS: 900 GeV 



(d) CMS: 900 GeV 



Covariance Matrix 



| 0.4149 -0.0008255 0.0006588 

1 | -0.0008255 2.737e-06 -3.047e-06 

2 | 0.0006588 -3.047e-06 4.026e-06 



X 2 = 26.364 
NDF (Degrees of freedom) = 20 
Z 2 ed = Z 2 /NDF = 1.318187 
Correlation = -0.9177 



0.066 0.068 0.07 0.072 0.074 

T [GeV] 



Covariance Malrix 



1.138; 
1.137; 

i.i3e; 
1.135; 
1.134; 
1.133; 
1.132; 
1.131; 

1.13: 



0.6096 -0.001058 0.0006346 
-0.001056 2.707e-06 -2.177e-06 
0.0006346 -2.177e-06 2.093e-06 



X 2 = 73.282 
NDF (Degrees of freedom) = 30 
X 2 ed = x 2 /NDF = 2.442722 
Correlation = -0.9147 



0.07 0.072 0.074 0.076 0.078 

T [GeV] 



(e) CMS: 2.36 TeV 



(f) CMS: 7 TeV 



Figure 3: Analysis of the correlation between the parameters T and q. x 2 is 
shown as a function of T and q and the cross show the position of the best 
fit values for the parameters. 
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ALICE: 0.9TeVx 10 4 
RHIC: 0.2TeVx 10~ 10 
ATLAS: 0.9TeVx 10~ 8 
CMS: 0.9TeVx 10~ 6 
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CMS: 7TeVx 10° 
T = 70 MeV 
T = 80 MeV 
T = 90 MeV 



10 



21 I I I I I 



I I I I I I I I I 



I I I I I I I 



10 



12 



14 



i i i I i i i 



1 p T (GeV/c) 



Figure 4: Typical results for fitting of Eq. [6] to experimental data with only 
q as a free parameter. The calculations are performed for three different 
temperatures. 



ranging from 0.2 TeV up to 7 TeV with a fixed temperature, T, and using 
only q as free parameter. Also, the new procedure allows us to study q as a 
function of \/i for different values of T. 

In Fig. Owe show q vs a/s as obtained for different values of T. We observe 
that q increases monotonically with y/s, the shape being approximately de- 
scribed by a sigmoidal function of the temperature. A sigmoidal behaviour 
was already conjectured in Ref. [3], but more experimental information is 
needed before drawing any canclusion, mainly at low energies. 

From the results shown in Fig. 2] we see that the best fittings are obtained 
for T 80 MeV, and in Fig. |5]we observe that for T = 80 MeV the entropic 
paramenter is approximately constant from yfs = lGeV to 7 GeV, with 
q « 1.12. 

The results obtained here are in agreement with an analysis performed 
in [5] , where the ]>r-distributions for different particles produced in p + p 
collisions at center-of-mass energy of 0.9 TeV were studied, resulting in T 
and q constant for all particles with T w 75 MeV and q ~ 1.15. There are 
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Figure 5: The values for q corresponding to the best fit with fixed T, for 
different values of temperature. Lines represent the best fitted curves for a 
sigmoidal function to the data. 

also indications that these results hold for A + A collisions [IB] . 

4 Conclusions 

In this work we present a systematic analysis of pr-distributions observed in 
p + p experiments. It is shown that the experimental data gives support to 
the hypothesis of a limiting effective temperature, T, and a limiting entropic 
factor, q. 

In the analysis, we show that T and q are correlated parameters in the 
fitting procedure of the theoretical transverse momentum distribution to the 
experimental data. From this correlation, it results that the critical temper- 
ature is T Q =(192±15) MeV, a value in good agreement with lattice-QCD 
predictions. 

The study presented here gives evidences for a limiting effective tem- 
perature for the hadronic system formed in hadron-hadron collisions, with 
T ~ 80 MeV, and also of a limiting entropic factor with q m 1.12. These 
results are in good agreement with values found in Literature. 



8 



5 Acknowledgements 

The authors thanks Prof. Dr Otaviano Helene (Universidade de Sao Paulo) 
for interesting discussions about the statistical analysis of the experimental 
data presented here. This work received support from the Brazilian agency, 
CNPq, under grants 305639/2010-2 (A.D.) and 136988/2009-1 (I.S.). 

References 

[1] R. Hagedorn, Suppl. Nuovo Cimento 3 (1965) 147. 

[2] I. Bediaga, E.M.F. Curado and J.M. de Miranda, Physica A: Statistical 
Mechanics and its Applications 286 (2000) 156 - 163. 

[3] C. Beck, Physica A: Statistical Mechanics and its Applications 286 
(2000) 164 - 180. 

[4] A. Deppman, Physica A 391 (2012) 6380-6385. 

[5] CMS Collaboration, Phys. Rev. Lett. 105 (2010) 022002. 

[6] ALICE Collaboration, Phys. Lett. B693 (2010) 53-68. 

[7] ATLAS Collaboration, Phys. Lett. B688 (2010) 21-42. 

[8] STAR Collaboration, Phys. Rev. Lett. 91 (2003) 172302. 

[9] J. Cleymans and D. Worku, J. Phys. G: Nucl. Part. Phys. 39 (2012) 
025006. 

[10] G. Wilk, Braz. J. Phys. 37 (2007) 714-716. 

[11] G. Wilk and Z. Wlodarczyk, Phys. Rev. C 79 (2009) 054903. 

[12] G. Wilk and Z. Wlodarczyk, Cent. Eur. J. Phys. 10 (2012) 568-575. 

[13] J.L. Rodgers and W. Alan Nicewander, The Am. Stat. 42 (1988) 59-66. 

[14] S. Ejiri, Nucl. Phys. B - Proc. Supp. 94 (2001) 19 - 26. 

[15] K. Miura, T.Z. Nakano, A. Ohnishi and N. Kawamoto, Phys. Rev. D 80 
(2009) 074034. 

[16] I. Sena and A. Deppman, arXiv:1208.2952vl (2012). 



9 



